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Abstract 

We study the statistical properties of spatially averaged global injected power fluctuations for 
^N I Taylor-Couette flow of a worm-like micellar gel formed by surfactant CTAT. At sufficiently high 

Weissenberg numbers (Wi) the shear rate and hence the injected power p{t) at a constant applied 
stress shows large irregular fluctuations in time. The nature of the probability distribution function 
CSJ . (pdf ) of p{t) and the power-law decay of its power spectrum are very similar to that observed in 

recent studies of elastic turbulence for polymer solutions. Remarkably, these non-Gaussian pdfs 
^^ . can be well described by an universal large deviation functional form given by the Generalized 

"^ ■ Gumbel (GG) distribution observed in the context of spatially averaged global measures in diverse 

C^ ■ classes of highly correlated systems. We show by in-situ rheology and polarized light scattering 

experiments that in the elastic turbulent regime the flow is spatially smooth but random in time, 
'^ ■ in agreement with a recent hypothesis for elastic turbulence. 

o . 

Q , PACS numbers: 

> 

in 
i> 
in 



o 



X 



Hydrodynamic instability like turbulence is an inertia driven phenomenon occurring at 
high Reynolds numbers (Re). However, for viscoelastic fluids of long polymers, similar 
instabilities have been observed even at a very low Re because of elastic hoop stresses 
generated by the stretching of the polymers in the curvilinear flow field, indicated by high 
Weissenberg number (Wi) (l| . Wi defined by the ratio of time scale set by the relaxation 
time (tr) of the fluid to that set by the strain rate (7), Wi = tr 7= Ni/a, with A^i is the 
first normal stress difference at a shear stress a, plays the same role for elastic turbulence 
that Re plays in the case of inertial turbulence. The elastic turbulence [l| is a flow instability 
that occurs at practically zero Re (Re<<l) and high Wi (Wi>>l). 

In past few years the statistics and scaling properties of the injected power fluctuations 
have been studied in detail both experimentally and numerically for inertial turbulence |2|-l7| 
in von Karman swirling flows. It has been shown that the probability distribution function 
(pdf ) of spatially averaged global injected power shows deviation from the Gaussian nature, 
skewness being on the side of smaller value of the injected power. Inertial turbulence also 
shows spatially smooth but random in time flow behaviour like elastic turbulence, below the 
Kolmogorov dissipation scale. Recently, similar results have emerged for power fluctuations 
in polymers in the context of elastic turbulence in von Kerman flow [8|. Although, the sign 
of the asymmetry of the power fluctuations appears to be same for both the inertial and 
the elastic turbulence, the magnitude of the normalized third moment or skewness vary in 
different studies. For some cases like in |6| the pdfs are very skewed (skewness ~ -1) but 
in many other cases [71, \^ this value is ~ -0.2. It has been argued that averaging over 
many independent large scale structures decreases the skewness of the global quantities. 
Very recently, the flow behaviour of worm-like micellar system has been shown to have 
signatures of elastic turbulence J9|. To our knowledge, there is no similar study of statistics 
and universality properties of spatially averaged global injected power in the case of unstable 
flow at high Wi and very low Re for worm-like micellar gels. These studies are important for 
worm-like micellar systems, since worm-like micelles have additional relaxation mechanisms 



due to scission- reattachment dynamics compared to polymers ^0[. This motivated us to 
study the spatially averaged global injected power fluctuations for Taylor- Couette flow of a 
shear thinning worm-like micellar solution which shows elastic turbulence. At sufficiently 
high Wi, the time-series of injected power p{t) = a 7 (t)K (where a is the apphed stress, 
7 (t) is the resulting shear rate and Vg is the volume of the sample) at a fixed shear stress 
shows large irregular fiuctuations. The pdfs of p{t) deviate from the Gaussian nature with 
skewness toward low power side with respect to (pit)) and the values of skewness obtained 
are similar to that in Newtonian fiuids [7| and polymer fiows pj despite the differences 
in the flow proflles in von Karman and Taylor-Couette geometries. Moreover, the three 
systems discussed above, namely, a Newtonian fluid [7|, polymer solutions |8| and the present 
worm-like micellar gel follow different routes to turbulence from the laminar flow because of 
their completely different dynamics and relaxation mechanisms which make their statistical 
similarities in the turbulent regime even more non-trivial. We show for the flrst time, 
that the non-Gaussian pdf of the injected power p{t) for elastic turbulence in worm-like 



micellar system can be described by an universal large deviation functional form given by the 
Generalized- Gumb el (GG) distribution, observed for fluctuations of injected power in inertial 
turbulent flows p, |ll|, [l2| and many other correlated systems 13l42l|. Physically, when the 



correlation length becomes comparable to the system size, the system can not be divided into 
a large number of independent domains for any global or spatially averaged measurements, 
hence, central limit theorem no longer holds and deviations from the Gaussian statistics is 
observed. To see the existence of this correlation directly, we performed in-situ polarized 
light scattering and rheology experiments as described in detail in |i22|. In the turbulent 
regime, the scattered light signals obtained from the gap of the transparent Couette cell 
show random fluctuations. A calculation of spatial and temporal correlation functions of 
the scattered intensity fluctuations reveals that although the time correlation is very short 
ranged, the system has very long ranged spatial correlation, practically limited by the size 
of the Couette cell itself. By this direct method we show the non-trivial space-time flow 
behaviour satisfying the criterion for the elastic turbulence in micellar system. To test the 
statistical robustness of the results, we repeated all the measurement in another Couette 
geometry with completely different aspect ratio. Both sets of results are in very good 
agreement, with each other. 

We used CTAT 2 wt% + 100 mM NaCl in water which forms a worm-like micellar gel 



for our experiments. The preparation of the sample is discussed in detail in [23|. All the 
experiments were carried out on a MCR300 stress controlled rheometer (Anton PAAR, Ger- 
many) at a temperature of 26 "C. For imaging experiments, we used home made transparent 
Couette geometry having inner cylinder diameter of 23 mm and height 40 mm and with a 
gap of 1 mm (Couette A). The outer cylinder is made of transparent glass and are partially 
enclosed by a water circulation chamber to control the temperature (the exposed portions 
of the glass cell were used for imaging). A thin sheet of polarized laser beam, formed by 
placing a cylindrical lens and polarizer arrangement in front of a randomly polarized He-Ne 
laser (10 mW) was sent along the vorticity (V x v) direction to illuminate the entire gap in 
(V X V, Vf ) plane and imaging was done perpendicularly to the illuminated plane with a wide 
angle lens and a CCD camera (Lumenera 0.75C, 640 x 480 pixels). By this arrangement 
the polarization of the incident laser beam can be tuned continuously in (v, Vf ) plane. All 
the optical components were mounted on XYZ stages for precision controls. For good sta- 
tistical accuracy, we analyze the shear-rate/stress time series with almost 50000 data points 
(data points during the start-up transients were ignored) sampled at a frequency of 5 Hz. 
The experimental tolerence in the value of applied shear-stress is ~ 10~^% and for applied 
shear-rate (obtained through a feed-back mechanism in our stress controlled instrument) is 
-1%. 

The flow curve for CTAT 2wt% + lOOmM NaCl is shown in Fig. la. The system shows 
Maxwellian behaviour [10| with a single relaxation time r^j ~ 1 s (inset to Fig. la). In the 
non-linear regime, in situ imaging experiments do not show shear banding, but reveal Taylor 
vortex flow right from the onset of the shear-thinning plateau. Also, at 7~ 10, there is a 
sudden jump in the stress value, indicating the enhanced flow resistance because of the 
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FIG. 1: (a) Shear stress (a) as a function of shear rate (7) for CTAT 2 wt% + 100 niM NaCl. The 
inset shows the Storage (C) and loss (G") moduhi as a function of angular frequency to, the solid 
lines indicate the fit to the single relaxation Maxwell's model. Shear-rate as a funtion of time for 
different apllied shear stresses (b) a = 17 Pa, (c) a = 33 Pa. 



onset of elastic turbulence. Beyond this value of 7 the Taylor vortices get destabilized and 
move randomly manifesting in large fluctuations in shear rate at a constant applied shear 
stress and vice versa 
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We study the temporal dynamics of the system in the non- 
linear regime for different applied shear stress values of 17 Pa and 33 Pa. The shear-rate 
does not exhibit fluctuations at stress values less than 17 Pa. For stress value of 17 Pa 
which corresponds to onset of elastic turbulence, shear-rate shows intermittent behaviour 
indicating that the system becomes unstable and shows huge fluctuations as shown in Fig. lb. 
Fig.lc indicates that the shear-rate fluctuates randomly with time for a stress value of 33 
Pa. This is the elastic turbulence regime and all the further statistical analysis are done in 
this regime. The dynamics of shear-rate exactly correlates with the dynamics of the Taylor 
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vortices as observed in 

The pdfs n(p) of injected power (p) in the turbulent regime for different applied shear 
stress values a = 34 Pa, 35 Pa and 36 Pa are shown in Fig.2a. From Fig. 2a we see that with 
increasing a, the peak of the corresponding pdf shifts toward higher value of the injected 
power. In Fig. 2b the same pdfs are plotted in terms of reduced value of injected power 
y = ■^, where, AP = p{t) — {p{t)). Here, {p{t)) and F are respectively the mean and the 
standard deviation of p{t). It can be seen that all the curves collapse into a master curve 
which deviate from the Gaussian function with negative skewness. The dotted line is a fit 
to Gaussian function to show the non- Gaussian nature of the master pdf. We now show 
that the non-Gaussian pdf of injected power pjt) shown in Fig.2b is well represented by 



Generalized Gumbel distribution (GG) given by [25 



U{y) = ir(a)e'^[-^('^)(^+^('^))-^"'*^^^'''"l (1) 

The above equation has a single parameter 'a' which is again fixed by the skewness 
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FIG. 2: Pdfs n(p) and n(y) of injected power for different applied stress values in the elastic 
turbulent regime, cr = 34 Pa, 35 Pa and 36 Pa in terms of (a) raw and (b) reduced (defined in the 
text) injected power. The solid line in (b) indicates the simulated curve using GG distribution for 
a = 39 and the dotted line is the fit to Gaussian function, (c) Pdf n(y) of spike width for y >0 
(hollow squares) and y <0 (solid circles) for an applied stress (a) of 34 Pa. The solid lines indicate 
the exponential fits (n(AT) = Iloe^ ' '^) to the data. For clarity pdf for y <0 have been shifted 
by a multiplicative factor of 10. (d) PSD of injected power fluctuations for applied stress values 
of 32 Pa, 33 Pa, 34 Pa, 35 Pa and 36 Pa (for Couette A) and also for 60 and 62 Pa (for Couette 
B, described in the text) The solid straight line indicates a power-law exponent of -4. The curves 
for Couette B are shifted by a multiplicative factor of 100 for clarity, (e) PSD of injected power in 
constant shear-rate mode for 12s~^, 13,s~^, 14s~^, 15s^^ and 16s^^ (Couette A). The solid straight 
line indicates a power-law exponent of -2. 



of the distribution as (y^) = l/y/a. All other parameters are related to 'a' as follows, 

b{a) = — ^^2 ; z{a) = i(/na ^^;7^^^) and K{a) = M^ where T{a) is the Gamma function. 

The non-Gaussian pdf with negative skewness matches perfectly to the simulated curve 
given by the GG distribution (Eq.l) as shown by the solid lines in Fig.2b over ~ 3 orders 
of magnitude with the value of the parameter a = 39, estimated from the third moment of 
the experimentally measured reduced injected power {y^). Fig.2c shows the distribution of 
temporal spike width (AT) (persistence time of a single fluctuation event) for the reduced 
injected power y above and below its mean value {{y) = 0) for an applied stress of 34 Pa. 
For both the cases n(AT) shows exponential tails over two orders of magnitude as shown 
by the linear fits in log-linear scale, indicating the statistical independence of the large time 



events, although for short times they remain strongly correlated. The time constants for 
the exponential decays are r = 2.2 s for y > and 1.4 s for y < 0. The mean spike width 
(AT) = 2.8 s for y > and 2.6 s for y < 0. Since, the skewness of n(j/) < (Fig. 2b), the 
higher values of (AT) and r for y > compared to those for y < implies that the large 
intermittent negative fluctuations have shorter life-time. The skewness of the distribution on 
the negative side has also been seen in polymer solutions 28| which has been postulated 
to be due to accumulation of the elastic stress near the static boundary. This observation 
is counterintuitive to the fact that that the shear-rate is higher near the inner cylinder 
and hence would result in larger stresses near the inner cylinder. More work is needed to 
understand our results. 

The power spectral density (PSD) of p(t) (i.e. 7 (t)) at different applied stress values 
in the turbulent regime are shown in Fig. 2d, clearly displaying a power- law decay with an 
exponent /3 = -4 over five orders of magnitude and a smeared peak at ~0.15 Hz. The value 
oi (3 satisfy the Fouxon [27| criterion (|/3| >3) for the elastic turbulence. To see the statistical 
robustness of the exponent /3, we repeated all the experiments in another Couette geometry 
having different aspect ratio: inner cylinder diameter of 32 mm and height 16.5 mm and 
a gap of 2 mm (Couette B). Here, the onset of elastic turbulence occurs at much higher 
value of applied stress (~ 50 Pa). The difference in the onset of elastic turbulence can come 
from the difference in curvatures of the two geometries. The value of /3 obtained for the 
PSD evaluated at the applied stress values of 60 and 62 Pa is again -4, as shown in Fig.2d. 
The PSD oi p{t) (i.e. cr(t)) at constant shear-rate for Couette A is shown in Fig.2e where 
again a power law with a slope of -2 is observed which does not satisfy the Fouxon criterion. 
We do not understand the origin of different power laws seen in controlled shear stress and 
controlled shear rate experiments. 
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FIG. 3: Space-time behaviour of (a) normalized spatial correlation function Cs{r) and (b) nor- 
malized temporal correlation function CT{t). (c) Time averaged spatial Fourier-amplitude I(k) of 
scattered intensity as a function of wave vector k (=1/1), the solid line is a power-law fit with 
exponent -1.8 to the data. The applied stress o" = 33 Pa in all cases. 



To look for any long range correlations in the system in the turbulent regime, we did 



in-situ polarized hght scattering studies as described in detail in [22|. We image the gap 



of the Couette cell A as a function of time in (V x v, Vv) plane in the turbulent regime 
with a CCD camera at a frame rate of 1 Hz. This frame rate is sufficient to capture the 



correlated dynamics in the system, since the injected power p{t) has a correlation time of 
~ 3 s. The image of the entire gap (~ 40x1 mm), is represented by a slice of 370x11 
pixels and we analyze 250 such slices in the turbulent regime. We study the space-time 
distribution of scattered intensity /(r, t), by studying the intensity along a vertical cut which 
covers the entire length of the cell along the vorticity direction (370x1 pixels) chosen near 
the center of the gap of the Couette cell and then stacking them at different times as the 
Space-Time Plot (STP). The STP of the scattered intensity distribution in the turbulent 
regime is quite random (data not shown). The normalized spatial correlation functions 
Csi'f') = {I{r')I{r + r')),^j evaluated at different times in the turbulent regime are shown 
in Fig. 3a which clearly reveals that the scattered intensity is correlated spatially over the 
entire length of the Couette geometry and practically limited by the size of the container, at 
every t. Similarly, we also calculate the normalized temporal correlation function Cxit) = 
{I{t')I{t + t'))j, for all the pixels along the cut. The correlation function Cxit) is extremely 
short lived, almost delta-correlated (up to our time resolution) as shown in Fig.3b. The 
temporal correlation function monotonously goes to the base-line over ~ 1 s as has been 
confirmed by a higher image grabbing rate of 3 Hz (not shown) . These observations clearly 
indicate that although the spatially averaged injected power fluctuates randomly in time, 
the spatial structures remain intact and they move randomly in time. This spatially smooth 
and temporarily random flow behaviour supports a recent hypothesis of elastic turbulence, 
not reported so far, in the context of micellar gels. Fig.3c shows the plot of time averaged 
spatial Fourier-transform of scattered intensity along the vertical cut mentioned above which 
show a power law decay (I(k) ~ k") with exponent a = —1.8 (similar to that obtained in 
[9|]) over a length scale spanning ~ 1 to 10 mm. 

In conclusion, we have studied the universality and scaling properties of spatially averaged 
global injected power p{t) in the context of elastic turbulence in a shear-thinning worm-like 
micellar gel. The non-Gaussian probability distribution functions (pdf) show a negative 
skewness very similar to the case of polymer solutions [8|]. The Power Spectral Density 
(PSD) oip{t) shows a power law behaviour with a decay exponent /3 = -4 (< —3), a signature 



of elastic turbulence [27j. But the decay exponent (3 is very different (~ —2) for constant 
shear-rate run which is not understood. We showed a universal scaling of the non-Gaussian 
pdf of p{t) in terms of Generalized Gumbel (GG) distribution, observed in many cases of 



spatially averaged global measures for a diverse class of highly correlated systems |13H21 
In-situ polarized light scattering studies clearly show a presence of spatially long ranged 
correlated dynamics in the turbulent regime although the temporal correlations are short 
ranged. We believe that the presence of this long ranged spatial correlation is responsible for 
the GG distribution of the spatially averaged global injected power. This kind of universality 



is similar to that observed for inertial turbulence p, [ll|, |l2| . The intermittent nature of the 
injected power below the mean value as compared to its value above the mean has been 
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postulated to be associated with the nature of the elastic stress field in the system 
We believe that our work will motivate further experimental and theoretical studies on the 
fascinating phenomenon of elastic turbulence in soft matter. 
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